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We describe, at the microscopic level, the dynamics of N interacting components 
where the probability is very small when N is large that a given component interact 
more than once, directly or indirectly, up to time i, with any other component. Due to 
this fact, we can consider, at the macroscopic level, the quadratic system of differential 
equations associated with the interaction and establish an explicit formula for the 
solution of this system. We moreover apply our results to some models of Econophysics. 

1 Introduction 

In their paper, Duffic-Sun (2012) (see also DufSe (2012)), the authors provide a 
mathematical fomidation for independent random matching of a large popula- 
tion. Here we develop an approach, inspired by Kac (1956), where the random 
matching is instead asymptotically independent. To do so, we start with a se- 
quence of dynamical sets of interacting components, one for each integer N. For 
these dynamical systems we can show that when N is large the probability is 
very small that a component has interacted more than once directly or indirectly 
up to time t with any other component. Thanks to this fundamental property, 
we can link the microscopic and macroscopic levels using results from the theory 
of continuous-time Markov chains. 

In section 2, we describe these dynamics with their symmetric interaction 
kernels. We consider interactions involving m components, for m > 2, and we 
suppose that the intensities of these dynamics have an adequate dependence on 
N. Our techniques enable us to obtain an explicit formula for the associated 
quadratic system of diffferential equations. We thereby extend the results first 
obtained in Duffie-Manso ( 2007) and pursued in Duffie-Giroux-Manso (2010). 
We note that our formula is valid for any interaction kernel and it is more explicit 
than the one obtained for the particular kernel considered in the latter article. 
It enables us, in particular, to obtain new results for models of Econophysics. 
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The statement and the proof that our formula solves the system of differential 
equations are done in section 3. In section 4, we present our applications to 
Econophysics. While in section 5, we come back to the fundamental property 
of the system's dynamics and obtain several intermediary results leading to its 
solution. 

2 The dynamics. 

We suppose that all components take their values in a measurable space, {E, £), 
(one can think of (R'^, and their interactions are given by a symmetric 

probability kernel Q on the product space (^E"^,£®"^ ) for m > 2. That is: the 
function Q{x\,X2, ■■■Xm',C\ x • • ■ x Cm) is measurable in {x\,X2, ■■■Xm)', is a 
probability measure in (Ci x • • • x Cm); and satisfies Q(xi, X2, ...x„i; Ci x ■ • • x 
1 •^cr(2)) ■■■^(7{m)'i ^ct(i) ^ ^f(2) X • • • X C^f^m)) fo'^ any permutation 

cr of {1, 2, to}. 

For each integer N ^ we consider an interacting set of N components which 
interact by groups of to according to the kernel Q. The interactions occur at each 
jump of a Poisson process with intensity A^. Groups are undistinguishable so 

each group has a probability of (^) ^ of being involved in a given interaction. 

We show, in section 5, the fundamental property that enables us to obtain 
the system's solution. In the next section, we describe the formula and show 
that it is indeed the system's solution. 

3 The solution of the macroscopic system. 

The kernel Q allows us to describe the macroscopic evolution of the system with 
an associated system of quadratic differential equations via the evolution of the 
law of a component. This probability law, denoted evolves with time and is 
in fact the solution of the Cauchy problem: 

=l^t - ; Mo = M 

where 



yLt°'"(C) = J ii{dxi)ii{dx2)-.n{dxm)Q{xi,X2,.-Xm;C X E"^ ^) ioi C & £ . 

The probability law /x°"'is the law of a component after the interaction of 
TO i.i.d. components with law /i. We can think of it as the law at the root of the 
TO-ary tree with only one interaction. We will look at all the trees representing 
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the interaction history of a component up to time t. So for a tree, A, with more 
than one interaction, we divide the tree in m subtrees at that last interaction 
and continue recursively up to time to define . (Please see figure 1 

of section 5 for an example of an interaction tree.) Let A„ be the set of all 
trees with n interactions (a.k.a. nodes), each node producing m branches. If 
■^n G A.„, then ^°mAn denotes the law obtained by iteration of /x""* through the 
successsive node of the tree when we place the law /i on each leaf of An- 

Now we will show that our Cauchy problem has a unique solution which can 
be expressed, by conditioning on the number of interactions up to time t, and 
then by the component's history. Such conditionings give us 

n-1 

where = ]^ + is the number of trees with n nodes, taking 

fc=i 

into account their branching orders; and Pn{t) is the probability of having n 
branchings up to time t. 

Finally, in order to show that the countable convex sum (1) gives the solution, 
we need the following lemma which will be proved in section 5. 

Lemma 1 = (*rff„, e"*(l - e-(™-i)*)". 

Theorem 2 The convex combination, 

n>0 *rn{n) ^^^^^ 

is the solution of the Cauchy problem 

= -Mt ;mo = 

Proof. Since the countable convex sum (*) is normally summable, we can 
differentiate /Uj term by term to obtain: 

-.. + e--'i:..(l-e-"-'T- E 

n>l ^ / \ / A„eAn 

Thus we need to show that: 

(C) ^ E n{l-e-^"^-'^r-' . _ ' -i^, E /^""^--(C) (2) 

n>0 ^ ' A„+ieA„+i 

Starting with the definition (on page 3), we have that the LHS of (2) is equal 
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to 



r™ Vi>0 ^ A,j6A,j / 



...Q(xi,...,x„;Cxi;™-i) 



which is equal to 



. _^ (m- l)"«i!...«„J 

which in turn is equal to 

i V^o (m-l)"n! y 

where 

F(ii,...,i„,n,/i, Aij,...., Aj^,(5,C) = ... 



And this last expression is a decomposition of the trees A^+i € A„+i appearing 
in the RHS of (2) in m subtrees after the first node (taking the branching order 
into account). The two expressions are therefore equal and this proves the 
theorem. ■ 

Remark 3 This sectAon brings a simplification to the special case studied in 
section 3 of Duffie-Giroux-Manso (2010) where the term a(^^_i^(^n-i)+i can now 
be given explicitly as (m-i)^-i{n-i)\ ■ ^^^^ suggests that some of the results of 
Duffie-Malamud-Manso (2009) can be extended to the case where information 
exchanges involve m agents. 
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Remark 4 We call the law = ^ e-*(l - e-("^-i)*)" ^ /x"-^" 

an extended Wild sum [11] and note that the convex combination we obtain for 
the case m = 2 is indeed the Wild sum, ii^ = ^ e^*(l — e~*)"^ /x""*"*", 

now well-known in the statistical physics of gases since the work of Kac (1956) 
[6J. 

4 An application in Econophysics: kinetic mod- 
els with random perturbations 



In Ferland-Giroux (1991) the authors study a class of kinetic equations of Kac's 
type and they show, for binary collisions, a convergence to the invariant law at 
an exponential rate. Bassctti et al (2011) show a similar result using different 
methods. The convergence in Ferland-Giroux (1991) is obtained along a set of 
convenient test functions with a telescoping technique due to Trotter (1959) and 
with the use of a version of Wild sums obtained from judicious conditioning. 
There is therefore the possibility of extending the results of Bassetti-Ladelli- 
Toscani (2011) for m > 2. We do this extension here only in the context of 
Ferland-Giroux (1991). For a review of Econophysics, in the case m = 2, one 
can read, among others, Diiring-Matthes-Toscani (2010) and Bassetti et al 
(2011) for more recent work. 

Let /i denote a probability law on M and let Smip) denote the law of HiXi + 
H2X2 + .■■+HmXm where the random variables {Xi] are independent and of law 
fi and the variables {Hi} are independant of each other and of the X variables. 
Under very general conditions, Rosier (1992) has shown that the transfoniation 
Sm has a fixed point. We will suppose here that the i?- variables have values 
in [0, 1] , that their mean is ^ and that they are not Bernoulli. We then have 
E[Hi + H2 + ... + Hm] = 1 and S'm(/i) has the same first moment as /i. If we 
suppose moreover that E[H^ + + ... + H^] < 1 then R6sler(1992) gives us 
the existence of a fixed point, denoted 7. This fixed point has a second moment 
as soon as does. Our goal is to establish the following result. 

Theorem 5 // we suppose that /i has a finite second m,oment then we have that 
the law = ^ e~*(l — e~('"~i)*)" (^_^^n„i ^ ^^Om-4„ converges to 7 at the 

n>0 A„eA„ 

exponential rate rj = 1 - E[H^ + H2 + ... -|- H^]. 

Our proof is an extension of Ferland-Giroux (1991) which treats the case m = 2. 
Proof. Let us first consider the tree with m leaves, denoted Ai. On each one 
of its leaves put independent random variables of law /i. Assume that these 
variables interact at a node to give HiXi + H2X2 + ... + HmX^- Let us call 
^o^Ai ^ Qj. jjjQj-g simply /Lt°"*, the law of this variable. In a similar fashion, we 
can consider 7°" (which is 7 since it is a fixed point). We will first consider the 



5 



differences | < , f > — < 7°",/ > | for each / G C^. One way to bound 
this difference is to use the telescoping technique of Trotter (1959) where we 
replace one by one (from the left say) the variables with law jj, by variables with 
law 7. We then obtain a sum of m terms of the form [ < 7°* o / > — < 

7°*=+! o > I which we will bound. In our particular models we can 

write down these expressions explicitly as 



\E [f{HiYi + H2Y2 + ... + HkYk + Hk+iXk+i + ... + H^X^)] - 

E [f{HiYi + H2Y2 + ... + HkYk + Hk+iYk+i + Hk+2Xk+2 + ... + HmX^)] \ 

where Yi : i = 1, ...k + 1 have law 7 and Xj, with j = k + 1, m, have law /i. 
Let Rk = HiYi + H2Y2 + ... + HkYk + Hk+2Xk+2 + ••• + -f^m^m- Then we have 

f{Rk + Hk+iXk+i) = f{Rk) + f'{Rk)Hk+iXk+i + f"{R*k) {Hk+iXk+if' 

and 

f{Rk + Hk+iYk+i) = f{Rk) + f'iRk)Hk+iYk+i + f"{RT) {Hk+iYk+if . 
Which in turn give us ■ 



<H°"^,f>-<J°-,f>\<cE 



,i=l 



We now need to iterate the process according to the different trees. Let An 
denote the set of leaves of the tree An . Then the contribution of a leaf u G An 
through its interactions down to the bottom of the tree is a product of the 
variables Hi with i = l,...,m. Let us denote this product by C„. For the tree 
An, the result of its interactions through the bottom of the tree will therefore go 
from CuXu, when all the variables put on leaves have law fx, to ^ C„F„ 

when all the variables have law 7. Applying the same techniques as above, 
namely a (longer) telescoping and a Taylor series expansion, we get 



J> - <l°"',f>\<cE 



u£A„ 



Let e„ = 



(m-l)i 



If we decompose An in its m subtrees 



from its first node we can apply a similar reasoning to Ferland-Giroux (1991) 
in order to obtain that e„ < cn""^ where a = and rj is such that 1 —rj = 



E 



E^f 



,i=l 



We now reformulate these assertions in propositions and and give their proof. 
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Lemma 6 We have e„ = (^) (l+ei + .... + e„_i). Therefore e„ < cn" ^ 



m— 1 * 



Proof. The decomposition of each tree An in m subtrees {AJj}™ ^ at the first 
node enables us to write e„ as a sum of m similar terms 



1 



(m - l)"n! 



.ueA' 



where is the set of leaves of the z*'' subtree. It suffices to treat the case 
1 = 1. Let us decompose A„, the set of trees with n nodes, by the number, 
fc, of nodes of the subtree A\. There arc ("j^^) (m — \)'^~'^~^{n — 1 — fc)! such 
trees. Indeed, since we need to take into account the order of appearance of 
these nodes wc have ("^^) choices for the appearances of A^'s nodes. Then we 
have (m — i)"-!-^ ways to divide the remaining nodes in the other m — 1 trees 
and finally, there are {n—l — k)\ choices for the nodes' appearances. Note that 
this number simplifies to (m — l)""^"*^. If we denote by A„,fc the subset 

of A„ formed by the trees An for which their subtree A\ has k nodes, we can 
then express /j as 



fi 



(m- 1) 



^ n — 1 

l)"n! ^ ^ 



fc=o A„eA„ 
21 n-1 



E 

.ueA'n 



E 



.u&A' 



= E[Hf\ 'f^ 1 ^ 
im - l)n (m - 1)^=^! ^^^^^ 

E\m] ^1 

= 7 — —r^ V efe , with eo = 1. 
(ra - l)n 

This proves the first assertion. A similar calculation to the proof of lemma 3 in 
Ferland-Giroux (1991) gives us the second result. ■ 



Proposition 7 For all f € we have | < /Xj, / > — < 7, / > | < ce 



Proof. Once again, it sufiices to follow the proof of theorem 3 in Ferland-Giroux 
(1991) with the extended Wild sum and replacing 1 — 77 by ■ 



5 The extended Wild sums. 



5.1 When graphs become trees 



In all our cases, we have an underlying market structure which is a Kac walk with 
interactions involving m agents. We add exponential times to obtain a marked 



7 



Poisson process whose marks are horizontal lines linking the agents participating 
in a given interaction. This enables us to describe the limit law of an agent, 
under an appropriate conditioning, as a countable convex combination on trees 
which is, as wc have shown in s(x;tion 3, the global solution of the associated 
differential equation on the space of probability laws. 

Here we explain how we came to that convex combination. We start our 
study by an analysis of the dynamics of the intrinsic structure of the large set of 
interacting agents when the number of agents increases. We assume that each 
interaction involves m agents, m > 2. More specifically, we consider a set of 

agents whose interactions happen at unexpected times so these interactions' 
occurrences follow a Poisson process. Since agents are interchangeable, each 

group has an equal probability of meeting of ( I . If we suppose the 

intensity of the meetings to be ^ then each agent has a meeting rate A which 
can be assumed to equal 1 under a time change. 

For A'' fixed and starting at time 0, we assign a vertical position to each 
agent. The down movement represents the passage of time, see figure 1 on page 
9. Each time a group of agents interacts, we draw a horizontal line between 
those agents and we draw a vertical line at each agent's position connecting to 
the horizontal line just drawn, so we see a random graph being formed. When 
we stop this graph at time t, we obtain the finite graph of all interactions that 
have taken place. Moreover, the history up to time i of a given agent, call it P, is 
described by the random sub-graph connecting all agents who have interacted 
directly or indirectly with P. A sample history of P's meetings/interactions 
may look like figure 1 below. 

The number of meetings is random but we can condition on it. The law of 
the finite graph is reversible since the meeting times are uniform on [0,t]. We 
want to show that a random graph representing the history of P can be replaced 
by a random tree as the number of agents, N, grows. If we look at figure 1, 
we see that the inclusion in the second meeting of one of the investors having 
participated in the first one, or the inclusion in the third meeting of an investor 
from the first or second one would create a cycle in our subgraph. As N grows 
though, the chance of meeting an investor previously encountered directly or 
indirectly tends to zero. 

To see this, let us consider the subgraph of P's history up to time t. Start- 
ing at time t, we pursue each one of the encountered vertical lines in P's 
history backward in time until wc reach the next horizontal line. If the in- 
clusion of the horizontal line in our graph does not create a cycle (i.e. no 
pair of investors were involved directly or indirectly in a previous meeting) 
we include the line, if not we remove it. Proceeding in this fashion up to 
time we get a tree with n branchings, say, which has the same law as the 
law of a tree obtained by a pure-birth process. Namely, the tree starting at 
P's vertical line at time t with intensity 1 and which at time has inten- 
sity (m — l)n -|- 1 and that same number of leaves. Between two branchings 
of this process a graph representing P's meeting history can have a random 
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Figure 1: 
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number of additional horizontal lines following a Poisson law of parameter at 

TO^^^™ 2^'^^^)(^m) ^' ^'^^ bound the expectation 

of these supplementary horizontal lines by a majorant which tends to as A'' 
increases. Indeed, since the mean number of redundant lines when there are 

n branchings up to time t is at most K ( f ~^ ^ ] ( ^ ] |)We 



2 J \ m ^ 

have that the mean number of redundant horizontal lines is bounded above by 

2''"^'^^ ) ( m ) j PJV,n(*)> where pjv,n(i) is the probability 

of having n branchings up to time t of the pure birth process with successive 
branching waiting times following exponential laws of parameter 

, N,, f N -{{m-l)n+l) \ f N 

A;v,„ = -((m-l)n + l)(^ m-1 )[m 



Since 



m V m-1 J \ '"T' 

_ ((m-l)n+l)r-»-f+^)) 

~ (N-i\ y"^) 

\m-l) 

< (m — l)n + 1 

then pNAi{t) is stochastically smaller than the law obtained with the intensities 
An = (™ — l)n + 1, which in turn are less than the intensities A„ = m{n + 1). 
Its transition kernel is then obtained by solving Kolmogorov's affine system of 
equations: 



mrvp^{n — 1) — m{n + l)p^{n) ; n > 1. 

Thus the latter intensities give us a geometric law Pf{n) = e"™*(l — e~™*)" . 
Since geometric laws have finite moments of all orders, the mean number of 
redundant horizontal lines is bounded above by a quantity converging to 0. 

For more details on Kolmogorov systems of equations for pure birth processes 
we refer the reader to Lefebvre (2006), for instance. 

Thus, after having specified the initial agents' states and their interaction 
kernels, we can approximate P' s law using the tree obtained from removing all 
redundant horizontal lines from its graph. We will use this fact in the next 
sub-section. 



dt 
dt 
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5.2 Limit countable convex combination 



We will now show that these random trees whose branching intensities depend 
on N can be approximated by trees with branching intensities independent of 
N. Taking into account that P's tree history is random with intensities de- 
pending on N, we could write P's law, denoted by /Xj'^, with complex formulae 
depending on N. Since our markets have a large number of investors, it is pre- 
ferrable instead to work with the limit of these laws. We note from (3) above 
that Ajv,n {{ni — l)n + 1). 

Let Pn{t) (— Pt{n)) be the solution of the afEne Kolmogorov system of equa- 
tions: 

dpt{n) 



= ((m - l)(n - 1) + l)pt{n - 1) - ((m - l)n + l)pt{n) ; n > 1. 
Recall fron the first section that ll^ = ^ Pn(0 #^\n) 1^°'"'^". 

n>0 A„eA„ 

Proposition 8 The sequence of laws yUj'^ converges to fj,^ as N increases. 

Proof. By Kurtz (1969) we have that PN.n{t) Pn{t) as N increases. But 
{Pnit))n>o is a probabihty law, so for e > 0, there exists n(e) such that 

p„{t) < e. Now let A''(e) be such that N > N{e) implies that \pN,n{t) - 

n>n(€) 

Pn{t)\ < ioi < n < n(e).We then have for C G f and N > N{e) 



n{e) 



*,N 



(C) - (C) I < ^ \pN,nit) - Pn it) I + 2e < 3e 



n=0 



since ^^(^-^ < 1 and {pN,n{t))n>o are probability laws. Our 

A„eA„ 

claim is proved. ■ 



Lemma 9 pt{n) = JrP„, e~\l - e"'™-!)*)' 



(m-l) 

Proof. We need to solve the afhne Kolmogorov system of equations (4). 
Proceeding by induction we have: 
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e-* 

t 

((m - l)(n - 1) + i)e-("(™-i)+i)* ^ e("('"-i)+i)%(n - l)ds 



To prove the lemma it suffices to note that ^m{n) = #m(n— l)((n— l)(m— 1) + 1) 
and that e("('"-i)+i)''e-^(l - e-(™-i)'*)"-i = e(™-i)^ (e^™"^)* - 1)""^ is the 
derivative of 7^-^ (e(™-i)« - 1)". ■ 

(m— l)n^ ' 

And this shows that the hmit law of P is indeed the extended Wild sum 
which we have shown to be the solution of the ODE associated to the interacting 
system . 

Acknowledgements. The first author would like to thank the Universite 
de Sherbrooke and its Faculte d'administration for their startup grant. 



dt 

dptjn) 
dt 
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